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Abstract 

We theoretically study multiphoton ionization through the triplet states of Mg by linearly po- 
larized (LP) and circularly polarized (CP) fs laser pulses. After the construction of the atomic 
basis using the frozen-core Hartree-Fock potential (FCHFP) as well as the model potential (MP) 
approaches for both singlet and triplet series which show rather good agreements with the existing 
data in terms of state energies and dipole matrix elements, we solve time-dependent Schrodinger 
equations with 3s3p ^Pi as an initial state, and calculate the total ionization yield and photoelec- 
tron energy spectra (PES). 
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I. INTRODUCTION 



During the last 30 years many theoretical and experimental investigations have been 
performed for Mg to obtain the atomic data and to understand its interaction with radiation 



through single-photon processes. 
of Mg were performed by Fischer 



^'he first extensive theoretical studies for the triplet states 
l| with a multi-configuration Hartree-Fock method which 
included correlations between the valence electrons, and by Victor and co-workers [2] with 
a semiempirical model potential which included core-polarization and dielectronic terms to 
calculate the oscillator strengths (OSs) for bound-bound transitions with ^'^S, ^'^P, and 
^'^D symmetries. Using a FCHFP with core-polarization and dielectronic terms Chang 
calculated the OSs between Ssnl ^'^L (L = S, P, D, etc.) states of Mg. Mendoza and Zeippen 
(4] studied photoionization from the excited triplet state 3s3p of Mg using a FCHFP 
with core-polarization and dielectronic terms in the close coupling approximation. Moccia 
and co-workers js] developed a nonempirical description of the core-polarization effects of 
Mg employing a basis set of modified Slater-type orbitals to study the transitions between 
the Ssnl ^'^L {L = S, P, D , and F) states of Mg. Luc-Koenig and co-workers j^l used an 
eigenchannel i?-matrix and multichannel quantum defect theory (MQDT) to investigate two- 
photon ionization of Mg atom. Lately, Fang and Chang [t] studied single-photon ionization 
from the excited singlet and triplet states of Mg below the Mg^+ threshold using an approach 
based on the 5-spline functions and Kim js] studied single-photon ionization from the 3s3p 
^'^P states with a i?-matrix method combined with MQDT. Most recently Fang and Chang 
has developed a i?-spline-based complex rotation method with spin-dependent interaction 
to calculate atomic photoionization of Mg with singlet-triplet mixing 9|. 

As for the multiphoton processes of Mg interacting with a laser pulse there are several 
experimental and theoretical works, all of which involve only singlet states: Kim and co- 
workers pi| studied single and double ionization of Mg by 10 ns Nd: YAG laser pulses at both 
532 and 1060 nm in the intensity range of 10^^ — 10^'^ W/cm^. Druten and co-workers [ll| 
measured PES associated with single and double ionization of Mg using 1 ps laser pulses in 
the wavelength of 580-595 nm and 10^^ — 10^'^ W/cm^ intensity range, respectively. Xenakis 
and co-workers |l^ investigated multiphoton ionization of Mg using 150 fs laser pulses at 
the wavelength of 400 nm for the peak laser intensities of up to 6 x 10^^ W/cm^. Gillen 

nn 

and co-workers [13|, IMj measured the ionization yield for single and double ionization of Mg 
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exposed to the 800 nm, 120 fs Ti:sapphire laser pulses for the peak intensities of 10^^ — 10^^ 
W/cm^, which was followed by the theoretical analysis 15 1. Liontos and his co-workers jiol 
investigated single and double ionization of Mg by Nd:YAG laser pulses with a ns duration 
for peak intensities up to 10^^ W/cm^. Zhang and Lambropoulos [l3] performed time- 
dependent calculations of Mg for the case in which ions are left in excited states. Recently 
we have studied the ionization yield and PES of Mg and clarified the origin of the subpeaks 



in the PES by the second and third harmonics of the fs Ti:sapphire laser pulse [18|. Note 
that all the previous studies have focused on multiphoton ionization from the singlet states 
of Mg. 

The purpose of this paper is to perform the theoretical study for the multiphoton ion- 
ization processes through the triplet states of Mg by LP and CP fs laser pulses. Specifically 
we choose 3s3p ^Pi as an initial state and perform time-dependent calculations after the 
construction of the atomic basis for both singlet and triplet series. This paper is orga- 
nized as follows. In Sees. HT] and UTTl we present the theoretical model: The time-dependent 
Schrodinger equation (TDSE), which describes the time- dependent interaction dynamics of 
the Mg atom with a laser pulse, is solved on the atomic basis states of Mg with two-active- 
valence electrons. Atomic units (a.u.) are used throughout this paper unless otherwise 
mentioned. In Sec. IIVI we present representative numerical results for the state energies and 
the J-independent and dependent OSs between the triplet states. Our results are compared 
with the existing data to confirm the accuracy of our atomic basis. Using those atomic 
basis states, we solve the TDSE to calculate the total ionization yield and PES from the 
3s3p '^Pi initial state of Mg by LP and CP fs laser pulses. Similar to the PES from the 
singlet ground state 3s^ ^5*1 of Mg [isl, the PES from the triplet 3s3p ^Pi state also exhibits 
subpeak structure. Finally, concluding remarks are given in Sec. |Vl 



II. ATOMIC BASIS STATES 



To start with, in order to study the interaction of the Mg atom with a laser pulse we 
have to construct the atomic basis of the Mg atom. The Mg atom is a two-valence-electron 
atom; it consists of a closed core (the nucleus and the ten inner-shell electrons ls^2s^2r)®) 

n 

and the two valence electrons. As it is already mentioned in the literature [19| there are 
several approaches to solve the Schrodinger equation for one- and two- valence-electron atoms 
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in a laser field. Since the general computational procedure has already been presented in 



Refs. 



2CH22j to construct the atomic basis states and the specific details about the atomic 



structure calculation of Mg have been reported in recent works [23|, |2J], we only briefly 
describe the method we employ. The field-free one-electron Hamiltonian of Mg"*", ha{r), is 
expressed as 



h(r)- ^ ^ + iii±ll + V M 



(1) 



where r represents the position vector of the valence electron, Z the core charge, / the orbital 
quantum number, and Ve//(r) the effective potential acting on the valence electron of Mg"*". 
Since the spin-orbit interaction is very weak for a light alkaline-earth-metal atom such as 



Me 
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261 ] ■ it might be safely neglected in the atomic Hamiltonian for our specific purpose. 



Similar to our recent study 2J] in which we have presented detailed comparisons between 
the frozen-core Hartree-Fock (FCHF) and MP calculations for the singlet states of Mg, we 
employ two different approaches in this paper to describe the effective potential, K//, in 
Eq. (H]). Namely (i) a FCHF potential and (ii) a MP. 



A. One-electron orbitals: Frozen-Core Hartree-Fock approach 

In the last years the most widely used method to describe the ionic core is the FCHF 
approach. In the FCHF approach the effective potential is given by 



VefAr) = yr{r)^Vr{r\ (2) 

where represents the FCHF potential and is the core-polarization potential which 
effectively accounts for the interaction between the closed core and the valence electrons 
2l| . Specifically we employ the following form for the core-polarization term: 



a. 



^AO = -^[l-exp-(^/'-')' 



(3) 



in which as = 0.491 is the static dipole polarizability of Mg^"*" j4| and r/ (/ = 0, 1, 2, ...) are 
the cutoff radii for the different orbital angular momenta: Tq = 1.241, ri = 1.383, = 1.250, 
rs = 1.300, and n = 1-100 H- 
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B. One-electron orbitals: Model potential approach 



Another simpler way to describe the ionic core is to use a MP, VJ^^ [2|, |24, l25|, |28| 
instead of the FCHFP, VJ^^. The advantage of the MP approach is that we can obtain the 
one-electron orbitals without self-consistent iterations, since the interactions of the valence 
electrons with the Mg^"*" core are replaced by pseudopotentials for each angular momentum. 
Thus the complexity of the problem is greatly reduced. That is, instead of the FCHFP, 
i.e., Vi^^{r), we employ the pseudopotential we have obtained in our previous work [2J] to 
describe the interaction of the valence electron with the Mg^"*" core: 

Vi^'^ir) = V/'{r) - ^exp {-ar^) + Biexp (-^r'), (4) 

where the values of the parameters introduced above, after the least-squares fitting, are 
A = 0.541, a = 0.561, Bq = 11.086, Bi = 5.206, 5/>2 = 0, /3o = 1-387, (3i = 1.002, 



s 



and 0i>2 = [2J]. We note that this form of is different from the one used in Refs. 



25^, 



In Sec. [TV] we will compare the results obtained by FCHFP, MP, and the 
experimental data. 

In either approach described above to obtain the one-electron orbitals, we employ a 
set of i?-spline functions to expand them. Thus solving the Schrodinger equation for the 
nonrelativistic one-electron Hamiltonian given in Eq. ([1]) is now reduced to an eigenvalue 
problem. 

C. Two-electron states 

Once the one-electron orbitals have been obtained using either the FCHFP or MP, we 
can construct two-electron states with the configuration interaction (CI) approach as we 
describe below: The field-free two-electron Hamiltonian, Ha{ri,r2), can be expressed as 

2 

Haivu rs) = + V{ru rs), (5) 

1=1 

where ha{ri) represents the one-electron Hamiltonian for the ith electron as shown in Eq. 
([1]), and l^(ri,r2) is a two-electron interaction operator, which includes the static Coulomb 
interaction l/|ri — and the effective dielectronic interaction potential Ti and r2 are 

the position vectors of the two valence electrons. By solving the two-electron Schrodinger 



equation for the Hamiltonian given in Eq. ([5]), the two-electron states are constructed with 



the CI approach 20N22|. For Mg, which is a hght alkahne-earth- metal atom, the LS coupling 
is known to give a good description and hence it is sufficient to label a two-electron state 
by the following set of quantum numbers: the principal, orbital, and spin quantum numbers 
for each electron, nj/jSj {i = 1,2), total orbital momentum L, total spin S, total angular 
momentum J, and its projection M on the quantization axis. After the CI procedure, two- 
electron states may be most generally labeled by the state energy and the quantum numbers 
iL,S,J,M). 

For singlet states (5 = 0), the above state labeling can be simplified to (L, M), since J is 
automatically equal to L. This is not the case, however, for the triplet states {S = 1) since 
J = L + S due to the presence of spin-orbit interaction. Physically, introduction of spin- 
orbit interactions influences the wave functions in two aspects: The dynamical (radial) part 
and the geometric (angular) part. As for the dynamical part we neglect its influence in this 
paper, since the spin-orbit interaction in the Mg atom is small j25|, anyway, as one can easily 
see from the very small fine structure splittings, and hence the radial wave function may be 
assumed to be J— independent as a lowest-order approximation. As for the geometric part, 
we can fully include it by introducing the additional quantum numbers, J and its projection 
M to specify the state. Thus it is necessary and sufficient that the triplet state is labeled 
hj{L,S,J,M). 

Now, once we have obtained the two-electron wave functions we are able to calculate 
the dipole matrix elements as well as OSs for both LP and CP fields. In the following 
two subsections we present two useful conversion relations between the J-dependent and 
J-independent dipole matrix elements and OSs, respectively. 

D. Calculation of the J-dependent dipole matrix elements 

By applying the well-known Wigner-Eckart theorem the following conversion relation 
exists between the J-dependent and J-independent dipole matrix elements if we define the 
initial and final states, i and /, by a set of quantum numbers 7j = {ui, L^, Si, Ji, Mi) and 
7/ = {uf, Lf, Sf, Jf, Mf), respectively: 

Dn.J.M.n,J„M, = (-l)^/-*^/+^/+^/+'^^+l+'^<55^,,^7(2 J, + 1)(2J; + 1) 
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/ Jf Si 
Ji Li 1 



DmL^ML^TifLfMLf^ (6) 



in which Dnij^M^nfjf,Mf ^"^^ L)n^L,ML nfLf,MLf represent the J-dependent and J-independent 
dipole matrix elements, respectively, q is associated with laser polarization, i.e., g = for 
LP and q = ±1 for right or left circular polarization (RCP or LCP), respectively. Recall 
that the allowed transitions take place between states accordingly to the dipole selection 
rules, which are generally written as J/ — Jj = 0, ±1 ( J/ — Jj = is forbidden if Ji = 0) 
and Mf — Mi = q {Mj = Mj = is forbidden if Jj — Ji = 0). In addition the following 
dipole selection rules are satisfied since L and S are good quantum numbers: Lf — Li = ±1, 
Mlj, — Ml^ = q and Sf — Si = 0, where Ml^^^^ represents the projection on the quantization 
axis of the orbital quantum momentum. 



E. Calculation of the J-dependent oscillator strengths 

Similarly, the OSs for multiplet transitions between two states, i and /, could be related 



to the J-dependent OSs 



29|: 



fM,nfJf) = {2L, + l){2Jf + l){ f{n,L,,nfLf), (7) 

' I Jf Lf] 

where f{niLi,nfLf) is the J-independent absorption OS, while f{niJi,nfJf) represents the 
J-dependent absorption OS. 



III. TIME-DEPENDENT SCHRODINGER EQUATION 

Having obtained the two-electron states constructed in a spherical box, we can now solve 
the TDSE. The TDSE for the two-electron atom interacting with a laser pulse reads 

2^^(ri, ra; t) = [Hair,, v^) + D{t)] ^'(ri, ra; t), (8) 

where \l/(ri,r2;t) are the total (two-electron) wave function at positions ri and r2 for each 
electron at time t, and Ha{ri,r2) is the field-free atomic Hamiltonian as shown in Eq. (JSj). 
The time-dependent interaction operator D{t) between the atom and the laser pulse is 
written in the velocity gauge as, 
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D(t) = -A(t)-(pi + p2), (9) 

where the dipole approximation has been employed, and pi and p2 are the momenta of the 
two electrons with A{t) being the vector potential given by 



A(t) = Ao/(t)cos(a;t). (10) 

Here Aq = A^qeg represents the amplitude of the vector potential and is the unit polar- 
ization vector of the laser pulse, expressed in the spherical coordinates, u and f{t) represent 
the photon energy and the temporal envelope of the laser field. In this paper we have as- 
sumed an envelope with a cosine-squared function, i.e., f{t) = cos^ (7rt/2r) where r is the 
full width at half maximum (FWHM) of the vector potential A(t). The integration time of 
Eq. dH]) is taken from — r to r. 

In order to solve Eq. ([8]), the time- dependent wave function, \l/(ri, t), is expanded on 
the atomic basis as a linear combination of two-electron states \l/(ri, En)'- 

vl>(ri, r2; t) = ^ Ce„jM(t)^(ri, r2; (11) 

nJM 

where CEr,jM{t) is the time-dependent coefficient for a state with an energy En, total two- 
electron angular momentum J, and its projection on the quantization axis M. Now, by 
replacing Eq. f|TT]) into Eq. ([8]) we obtain a set of first-order differential equations for the 
time-dependent coefficients CE„jM{t): 

i-T,CE„JM{t) = X! [En^nn'^JJ'^MM' — DnJMn'J'M'{t)\CE^,J'M'{t), (12) 
"''^ n',J',M' 

where DnjMn'j'M'if) represents the J-dependent dipole matrix element calculated in Sec. 
Ill Dl between two triplet states defined by the quantum numbers (nJM) and {n'J'M'). This 
means that we have neglected the spin-forbidden transitions between triplet and singlet 
states, which is reasonable for a light atom such as Mg. Specifically in what follows, we 
assume that the Mg atom is initially in the triplet state of the lowest electronic configuration, 
3s3p ^Pi{M = 0), i.e., 

\CE„jM{t = -r)!^ = Sn35jihio- (13) 
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The relevant energies of triplet states, averaged over the multiplet components, are presented 
in Fig. [TJ Note that 3s3p ^Pj is the lowest triplet state located at approximatively 2.71 eV 
from the ground state, 3s^ ^S. Our specific choice for the initial state results in the great 
simplification of the time-dependent problem to deal with, since the allowed transition paths 
by the LP or CP field become very simple as shown in Figs. [2]^a) and[2]^b), respectively. If 
we chose a different initial state, for instance, 3s4s ^Si, the transition paths for the LP field 
would be far more complicated than those shown in Fig. [2](a), and accordingly Eq. f lT2|) 
would become much more difficult to solve due to the enormous complexity of the transition 
paths. In contrast, this kind of complexity does not happen for the transitions between the 
singlet states jl8 |. 

Once we have obtained the time-dependent coefficients Ce„jm by solving Eq. f[T^ . the 
ionization yield Y and PES dP/dE can be calculated at the end of the pulse: 

Y=l- \CE^JM{t = +T)\\ (14) 

n,J,M{En<0) 

and 

dP 

= E \CE„jM{t = +r)\\ (15) 

"'^ E„=Ee J,M(En=Ee) 

where Eg represents the photoelectron energy. 



IV. NUMERICAL RESULTS 



Before solving the TDSE we must perform several checks regarding the accuracy of the 
atomic basis for the triplet states of Mg. Related to this, we have already obtained accurate 
atomic basis for the singlet states in our previous work [24] using FCHFP as well as MP 
approaches. The atomic basis states we need to solve the TDSE is constructed in a box size 
of 300 a.u. for the total angular momentum up to J = 9 with 1000 states for each total 
angular momentum. A number of 302 5-spline polynomials of order 9 with a sinelike knot 
grid is employed. To check the numerical convergence we have increased the box size up 
to 1000 a.u. together with an increased number of total angular momentum up to J = 14 
for each given intensity. It turned out that the basis states constructed in a box of 300 a.u. 
with the total angular momentum up to J = 9 with 800 states for each angular momentum 
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are sufficient to obtain a reasonable convergence in terms of the total ionization yield as well 
as PES. In Table I we present the two-electron angular configurations of type {rii li, n2 h) 
included in the construction of the two-electron wave functions. The principal quantum 
numbers are taken values in the range rii = (3 — 7) and 77-2 = (1 — 290) (with rii 7^ n2 if 
h = I2 ), respectively for each symmetry. The number of the two-electron configurations 
varies between 1100 and 1300 for the total angular momentum up to J = 9. 

To start with, we have compared the OSs for the triplet states obtained by the length 
and velocity gauges with the FCHFP approach, and confirmed that the agreement is quite 
good. This is a good indication that our wave functions are accurate. As for the MP 

approach, however, it is well known that the physically correct dipole matrix elements can 

ii 

be calculated only in the length gauge I30i] . since the Hamiltonian becomes nonlocal due to 



24l|). and we cannot perform a similar comparison 



the /-dependence of the MP (see Ref. 
between the two gauges. 

As a more direct comparison, we have calculated the state energies, OSs, and dipole 
matrix elements by both FCHFP and MP approaches and compared them with the existing 
theoretical and experimental data. In Table II we show the comparison of the calculated 
energies for the first ionization threshold and the first few triplet states 3snl "^L with the 
corresponding experimental values, where n = (3 — 6) or (4 — 7) for each total orbital 
momentum L = S,P,D and F. The energies (in units of eV) are taken with respect to 
the second ionization threshold Mg^"*" and the triplet states energies are averaged over the 
multiplet components. The theoretical data are taken from Ref. [jl and the experimental 
data are taken from the database of the National Institute of Standards and Technology 
(NIST) ^]. There is an overall good agreement between the calculated energies and the 
experimental values, and in addition our MP approach provides more accurate energies than 
our FCHFP approach. Of course, the accuracy of the energies do not guarantee the accuracy 
of the wave functions, and we must further check the accuracy of the wave function in terms 
of the J-independent and dependent OSs. 

Table HI presents the comparison of the J-independent OSs for single-photon transitions 
calculated by the FCHFP and MP approaches with other theoretical works Refs. |l|, 3], 0] 
and the experimental data taken from NIST. The OSs in the length gauge are shown for 
single-photon transitions among the first few triplet states: 3s4s — ?■ 3s(4 — 7)p ^P, 
3s3p 3p ^ 3s(4 - 7)s '^S, 3s3p 3s{3 - 6)d ^D, 3s3d -> 3s(3 - 6)p ^P, and 
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3s3d ^ 3s (4 - 7)/ ^F. From Table III it is clear that both FCHF and MP approaches 
provide an accurate atomic basis for the triplet states of Mg, and the overall agreement is 
quite well with other accurate calculations and the experimental data. There are, however, 
relatively large differences in the OSs for the 3sAs — 3s6p transition calculated by 
the FCHFP approach and 3s3d — t- 3sAp ^P transition calculated by the MP approach. 
Besides a small difference exists in the OS of the 3s3d ^D^ — )■ 3s4p ^P° transition calculated 
by both FCHF and MP approaches. This might be due to the very small energy difference 
between these two bound states of 0.016 eV. 

Finally, in Table IV we show the comparison of the J-dependent single-photon absorption 
OSs calculated by the FCHF and MP approaches with the experimental data taken from 
NIST. The calculated OSs are shown for the length gauge for the single-photon transitions 
among the first few triplet states: 3s4s ^5*1 — > 3s4p ^-Po,i,2 , 3s3p ^-Po,i,2 — ^ 3s4s '^5'i , 
3s3p 3Po,i,2 ^ 3s3d 2,3 , 3s3d 2,3 ^ 3s5j9 ^Pq ^ ^ , and 3s3d ^Di^2,z ^ 3s4/ ^^2,3,4 • 
Again, the overall agreement is quite good between our results and the experimental values. 
Therefore in what follows we present numerical TDSE results using the the atomic basis 
calculated by the FCHFP only. 

Having checked the accuracy of the atomic basis for the triplet states, we are now ready to 
perform the time integration of Eq. (fT2|) under various intensities for both LP and CP laser 
pulses. Recall that a number of 800 two-electron states for each total angular momentum 
up to J = 9 was used for the numerical integration of TDSE, thus leading to a total number 
of 7200 coupled differential equation to be solved. Please note that the typical size of the 
dipole matrices is about 800 x 800. The Runge-Kutta subroutines were used to perform the 
numerical integration of TDSE. As we have already mentioned, our initial state is 3s3p ^Pi 
(M = 0) and the photon energy is 2.7 eV which can be obtained from the second harmonic 
of a Ti:sapphire laser. Since the energy difference from 3s3p ^Pi to the ionization threshold 
is about 4.93 eV, at least two photons are needed for ionization. The intensity range we 
have considered for the numerical calculations is from 10^^ W/cm^ up to 10^^ W/cm^. The 
Keldysh parameter 7 is 1.1 at 10^^ W/cm^. 

The last check we should perform is that we may neglect the entire singlet states when we 
solve the TDSE for the triplet states. This check is particularly important, since our photon 
energy (2.7 eV) is resonant with the spin-forbidden 3s^ ^5*0 - 3s3p ^Pi transition. Because 
we cannot calculate the dipole matrix elements for spin-forbidden transitions within the 
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method we use, we have taken the experimental OS for the spin- forbidden ^5'o - 3s3p ^Pi 
transition, 2.38 x 10^^ a.u. from NIST, which is at least five and six orders of magnitude 
smaller than those for the (nearest) 3s3p — t- 3s3d and — > 3s3p ^P 
transitions, respectively. By phenomenologically including this spin-forbidden 3s^ ^Sq - 
3s3p ^Pi transition as shown in Fig. 121 we now solve two sets of TDSEs for the singlet 
and triplet series which are coupled through the resonant but very weak spin-forbidden 
3s^ —7- 3s3p ^P transition. After solving the two sets of TDSEs, we have ensured that, 
provided the 3s3p ^Pi initial state, the influence of the singlet states is extremely small as 
we expected, and we have safely neglected them in the following numerical calculations. 



A. Ionization yield 

The ionization yield is shown in Fig. Hl^a) with a log-log scale as a function of peak 
intensity for the LP (solid) and RCP (dashed) pulses. For the photon energy 2.7 eV we 
have chosen, both curves have a linear dependence on the peak intensity, up to 10^^ W/cm^, 
with a slope of 1.9, indicating that our results agree well with the prediction of lowest-order 
perturbation theory (LOPT). For peak intensities higher than 2 x 10^^ W/cm^, saturation 
starts to take place. Figure Hl^b) presents the ratio between the ionization yield by the CP 
and LP pulses, Ycp/Yip, as a function of peak intensity. For peak intensities up to 10^^ 
W/cm^, the ionization yield by the RCP pulse is about 0.83 times smaller than that by 
the LP pulse. Figures |l](a) and \Mjo) also suggest that ionization by the LP pulse is more 
efficient than the RCP pulse when ionization starts from the 3s3p ^Pi initial state. This 
result is somehow different from our previous time-dependent calculations for multiphoton 
ionization of Mg isl: It showed that, when less than four photons are needed for ionization, 
ionization from the singlet state with symmetry by the CP field starts to become more 
efficient than that by the LP field for a wide range of photon energy. In the LOPT regime 
the main reason that the ionization yield by the CP pulse is larger or smaller than by the 
LP pulse, for a non-resonant photon energy, is determined by the particular values of the 
total angular momentum and its projection on the quantization axis. 
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B. Photoelectron energy spectra 



In Fig. [5]we present representative results of the PES by the LP (sohd) and RCP (dashed) 
pulses at the peak intensity of 5 x 10^^ W/cm^. As it goes to the higher orders of above 
threshold ionization (ATI), the height of the ATI peaks by the LP pulse is more than one 
order of magnitude larger than that by the RCP pulse. Of course, this could be qualitatively 
understood that photoionization by the LP pulse has more chance to be near resonance with 
bound states than the CP pulse, and in addition there are more accessible continua for the 
LP pulse. It is interesting to note that subpeaks appear between the main ATI peaks, 
labeled as (b) and (c), for both LP and RCP pulses, and the height of the subpeaks is at 
least 5 orders of magnitude smaller than that of the main peaks. In addition there are small 
subpeaks, labeled as (a), on the right-side shoulders of the main peaks for both LP and 
RCP pulses. These results are reminiscent of the subpeaks studied in our recent paper for 



the singlet states of Mg in which multiphoton ionization of Mg from the singlet ground 
state has been theoretically studied. In that paper the origin of the subpeaks is clearly 
attributed to the bound states 3snp {n = 3,4,5...) which are far off-resonantly excited 
by the spectral wing of the pulse. In the next subsection we will identify the origin of the 
subpeaks in PES in a similar manner. 

Figures [6]^a)-l6]|^c) show the variation of the PES for three different pulse durations, (a) 
r = 80 fs, (b) 40 fs, and (c) 20 fs (FWHM). The photon energy and peak intensity are 
2.7 eV and 5 x 10^^ W/cm^, respectively. As the pulse duration decreases the ATI peaks 
are broadened and their heights are decreased. Besides, the subpeaks gradually disappear 
because of the broadening of the Fourier bandwidth of the shorter pulse. 

C. Origin of the subpeaks in the photoelectron energy spectra 

The method we have used to identify the origin of the subpeaks in the PES mentioned 

i' I 

in the previous subsection is quite similar to the one employed in our previous work [18] for 
singlet states of Mg: If the subpeaks arise from some photoionization processes involving 
four or five photons to leave the ionic core in some excited state, the height of the subpeaks 
with respect to the main peaks would be even much smaller than those in Fig. [S]at the peak 
intensity of 5 x 10^^ W/cm^, assuming the typical excitation/ionization efficiency with four 
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or five photons. The subpeaks cannot be attributed to some intensity-dependent effects, 
either; the ponderomotive shift is as small as 0.098 eV at peak intensity 5 x 10^^ W/cm^, 
and there are no triplet states coming into resonance during the pulse duration for both LP 
and CP pulses. Perhaps the subpeaks originate from the off-resonant excitations of some 
bound states, which, however, must be confirmed by the numerical calculations. Since we 
propagate the TDSE on the atomic basis, we can easily check this by solving the TDSE 
after the removal of the particular bound state under suspect, and comparing the PES with 



the original one with all states included 18 1. 



In Figs. [71^a)-[7I^c) we summarize the results for the PES calculated with the LP pulse at 
the peak intensity of 5 x 10^^ W/cm^. They are the results obtained after the removal of a 
particular bound state, namely (a) 3s3d ^Di, (b) SsAd ^Di, and (c) 3s5d ^Di, upon solving 
the TDSE, and compared with the result with the complete calculation of PES including all 
atomic triplet states of Mg. When the 3s3d ^Di state is removed [Fig. El^a)], the spike on the 
right-side shoulders of each main peak disappears. In addition the height of the main peaks 
is reduced since the 3s3d ^"Di state brings an important (but nonresonant) contribution to 
the ionization process. In this particular case the laser detuning is 0.53 eV with respect 
to the 3s3p ^Pi state. That is, the small spike, located at 1.02 eV, corresponds to the 
single-photon ionization process from the off-resonantly excited 3s3d ^Di state. Similarly, 
by removing the 3snd ^Di {n = A and 5) states different subpeaks labeled as (b) and (c) in 
Fig. [5] disappear, as can be seen in Figs. ^Ih) and[7l^c). This indicates that the physical 
origin of the subpeaks (a), (b), and (c) in Fig. [5]for the triplet states of Mg is quite similar to 
that we have found for the singlet states of Mg ISj. Briefly, off-resonant bound states such 
as 3snd ^Di {n = 3, 4, 5,...) are the origin of the subpeaks. Note that these states are located 
at 5.94, 6.71, and 7.06 eV, respectively, from the ground state 3s^ ^S, and accordingly the 
corresponding detunings are 0.53, 1.3, and 1.65 eV from the 3s3p ^Pi state since the photon 
energy is 2.7 eV. As for Fig. Wl^c) we note that, in addition to the 3s5d ^Di state, another 
not-identified state(s) might contribute to the subpeaks of interest labeled as (c). 

V. CONCLUSIONS 

In conclusion, we have theoretically studied multiphoton ionization of Mg from the triplet 
3s3p ^Pi state by linearly and circularly polarized fs pulses. For that purpose we have first 
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constructed the atomic basis with J— dependent dipole matrix elements for two active elec- 
trons, and then solved time-dependent Schrodinger equations with them. Since the spin-orbit 
interaction is rather weak for the Mg atom, J— dependent dipole matrix elements obtained 
by only taking into account the geometric (angular) part of the wave functions result in 
rather accurate values and compare well with the existing theoretical and experimental 
data. For the time-dependent calculations for multiphoton ionization from the triplet 3s3p 
^Pi state, the photon energy we have specifically chosen is 2.7 eV and corresponds to the 
3s3p ^Pi — > 3s^ ^Sq transition which is spin-forbidden and extremely weak. We have 
ensured that, even for the resonant photon energy, the singlet states do not influence the 
photoionization process. The ionization yields have been found to be larger for the linearly 
polarized pulse than for the circularly polarized pulse. Since the Mg atom has a rather 
rich level structure, the photoelectron energy spectra exhibits subpeaks in addition to the 
ordinary main ATI peaks. We have clarified the source of those subpeaks as ATI originating 
from some triplet bound states which are far off-resonantly excited by the spectral wing of 
the pulse. 
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TABLE I. Types of two-electron angular configurations used for the construction 
of two-electron wave functions. 
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TABLE II. Comparison of the energies for the first ionization threshold and the 
first few triplet states of Mg. The energies (in units of eV) are taken 
with respect to the second ionization threshold Mg^"*" and the triplet 
state energies are averaged over the multiplet components. 







FCHFP 


MP 


Theory [7] 


Exp(NIST) 


-E'Mg+ 


-15.000 


-15.042 




-15.035 




3ge 
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-15.281 
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TABLE III. Comparison of the J-independent single-photon oscillator strengths (in 
a.u. and length gauge) between the first few triplet states with ^5*^, 
^P°, ^D^, and ^F^ symmetry. Numbers in square brackets indicate 
powers of 10. 



3s4s '^S^ 


3s4p 


3s5j9 


3s6p ^P° 


3s7j9 
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1.320 
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6.963[-3] 
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MP 
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1.266[-1] 
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2.91[-2] 



20 



TABLE IV. Comparison of the J-dependent single-photon absorption oscillator 
strengths (in a.u. and length gauge) between the first few triplet states 
with ^Sj, ^Pj, ^Dj, and ^Fj symmetry. Numbers in square brackets 
indicate powers of 10. 



FCHFP MP Exp(NIST) 
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FIG. 1: (Color online) Energies of the triplet states of Mg. In order to show the relative positions 
with respect to the singlet ground state, 3s^ is also shown in this energy diagram. 
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FIG. 2: Allowed ionization paths in multiphoton ionization from an ^Pi (M = 0) initial state by 
(a) a linearly polarized field and (b) a right circularly polarized field. 
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FIG. 3: (Color online) Ionization scheme for Mg, at the photon energy of 2.7 eV, which includes 
the atomic basis of both singlet and triplet states resonantly coupled through the 3s^ ^S* — >■ 3s3p 
■^P transition. 
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FIG. 4: (Color online) (a) Ionization yield as a function of the peak intensity for linearly (solid) 
and right circularly polarized (dashed) laser pulses at the photon energy of 2.7 eV. The initial 
state of Mg is the triplet state 3s3p ^Pi (M = 0) and the laser pulse duration is 20 fs (FWHM). 
(b) Ratio of the ionization yield by the right circularly polarized pulse Ycp to that by the linearly 
polarized pulse Ylp- 
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FIG. 5: (Color online) Photoelectron energy spectra by the linearly (solid) and right circularly 
polarized (dashed) laser pulses at the photon energy of 2.7 eV. The initial state of Mg is the triplet 
state 3s3p ^Pi (M = 0). The pulse duration and peak intensity are 20 fs (FWHM) and 5 x 10^^ 
W/cm^, respectively. 
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FIG. 6: (Color online) Photoelectron energy spectra by the linearly laser pulses at the photon 
energy of 2.7 eV and the peak intensity 5 x 10^^ W/cm^ for tree different values of the pulse 
duration (a) 80 fs, (b) 40 fs, and (c) 20 fs (FWHM). The initial state of Mg is the triplet state 
3s3p ^Pi (M = 0). 
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FIG. 7: (Color online) Comparison of the photoelectron spectra by a linearly polarized laser pulse 
at the photon energy of 2.7 eV when the (a) 3s3d '^Di, (b) 3s4d ^Di, and (c) 3s5d ^Di bound states 
of Mg are removed from the atomic basis when solving the time-dependent Schrodinger equation. 
The pulse duration and peak intensity are 20 fs (FWHM) and 5 x 10^^ W/cm^, respectively. 
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